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ABSTRACT 



Aims. We investigate whether the so-called kappa distribution, often used to fit electron distributions detected in-situ in the solar wind, 
can describe electrons producing the hard X-ray emission in solar flares. 

Methods. Using Ramaty High Energy Solar Spectroscopic imager (RHESSI) flare data we fit spatially- and feature-integrated spectra, 
assuming a kappa distribution for the mean electron flux spectrum. 

Results. We show that a single kappa distribution generally cannot describe spatially integrated X-ray emission composed of both 
footpoint and coronal sources. In contrast, the kappa distribution is consistent with mean electron spectra producing hard X-ray 
emission in some coronal sources. 

Key words. Sumflares - Sun: X-rays, gamma rays - Methods: data analysis 



1. Introduction 

In solar flares, electrons are very efficiently acceler- 
ated, as shown by non-thermal X-ray and radio emission 
dTandberg-Hanssen & Em"shel ll988t iBenzll 19931 lKarhckvj[l997l 
iKlein & MacKinnonl2007l) . A widely used method to determine 
the electron distribution functions in solar flares is the forward- 
fitting of X-ray spectra in which the low-energy part usually 
is described by isothermal bremsstrahlung while for higher 
energies, electron pow er-law distributions with a low-energy 
cutoff are assumed (e.g. Hol man et all 2003). On the other hand, 
energetic electrons observed in-situ in space have a distribution 
with an enhanced number of particles in the high-energy tail. 
These electrons are commonly modelled by kappa distribution s 
jVandasll989HMace & Hellberdl995tlMaksimovic et alj 1997b . 
This leads us to the idea that the kappa distribution is a natural 
consequenc e of underl y ing ph y sical proce ss es. In a series of 
papers by lYoon et alj (12006b : iRhee et all d2006l) : iRvu et all 
d2007b it was shown analytically and numerically that the kappa 
distribution is the byproduct of beam-plasma interaction which 
first leads to the excitation of Langmuir waves followed by 
turbulent mode-coupling processes. Among the various wave- 
particle and wave-wave interaction processes, it was proposed 
that the "collisionality", defined via the number of particles 
per Debye sphere, plays a crucial role. Meanwhile, it appeared 
that the problem o f the kappa di s tribution has an even deeper 
physical meaning. iTsallisI (11988b : iLeubnerl (12002b have found 
that the kappa distribution is the thermodynamic equilibrium 
solution within so-called nonextensive thermodynamics with a 
generalised entropy definition. 

Although in the interpretation of the flare X-ray spectra the 
kappa distribution has not been considered up to now, it is evi- 
dent that it has some attractive properties: a smooth transition be- 
tween thermal and high-energy parts of the X-ray emission pro- 
ducing electrons and no need for any low-energy cutoff. Indeed, 
the recent statistical study by iKontar et alJ (2008) using a reg- 
ularised inversion indicates that if low-energy cutoffs exist in 



mean electron spectra, their energies are likely rather low (be- 
low 12 keV). Furthermore, taking into ac count velocity diffusion 
and the slowing-down rate of electrons. iGallowav et alj {2005) 
applied a steady-state solution of the Fokker-Planck equation to 
RHESSI spectra, thus consistently treating thermal and fast elec- 
trons. Finally, we note that the kappa distribution is also consid- 
ered to be connected to t he enhanced EUV line emission ob- 
served during solar flares (|Dz ifcakova et al. 200U). 

RHESSI dLin et al.l l2002) not only provides high resolution 
spatially integrated spectra but it is possible to obtain feature 
integrated spectra, i.e. spectra fr om individual sources, using 
imaging spectroscopy techniques (Battaglia & Benz 2006). 

The aim of this paper is to examine whether flare X-ray spa- 
tially integrated spectra or spectra of some particular type of 
source are consistent with a kappa distribution of the electrons 
producing the emission. We use the forward-fit approach and 
fit the spectra assuming a kappa distribution of electrons in the 
thin-target model. 

Properties of the kappa distribution and expressions for 
the corresponding bremsstrahlung thin-target emission are pre- 
sented in Section [2] The data analysis is described in Section [3] 
In Section|4]we show which types of flare spectra and of which 
sources are consistent with the kappa distribution of electrons. 
Our results are summarised in Section|5] 



2. Kappa distribution 

The kappa distribution f K (E) represents a generalisation of 
a Maxwellian distribution with an enhanced high-energy tail 
(ISummers & Thorne1ll991b . It can be written as 



ME) = A K - 
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where A K = r(K+l)/[F(K-0.5)(K-1.5) L5 ] so that J f K (E)dE = 1, 
T K is "temperature", k B is the Boltzmann constant, and E denotes 
electron kinetic energy. At large E and for low k, the distribution 
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is a power-law; in the limit k — > oo it approaches a Maxwellian 
distribution. 

2.1. Thin-target X-ray emission 

Assuming that superthermal electrons are described by the dis- 
tribution (Q]) and the mean electron density N K in optically 
thin plasma volume V, the thin-target bremsstrahlung photon 
spectrum 7(e) [ph otons cirT 2 s _1 per unit e] has the form 
dBrown et al.ll20Q3h 

7TVN f°° r 
I(e) = ~biFj A( £ M £ )2( e ' £ ) d£ 7 b = J npdV/V, (2) 

where Q(e, E) is the bremsstrahlung cross-section differential in 
photon energy e, v denotes velocity, and R is the distance to the 
emission source, e.g. 1 AU. Note that contrary to the widely used 
power-law forms for electron distributions producing solar flare 
bremsstrahlung emission, there is no need to introduce any (of- 
ten arbitrary) low-energy cutoff. Furthermore, the kappa distri- 
bution also describes the low-energy part of the X-ray emission 
that is usually fitted by a thermal component. 

3. Data analysis 

In order to explore a range of source types, we considered flares 
of a significant thermal component as well as flares with reported 
non-thermal behaviour at low energi es (~ 10 keV) - the so- 
called early impulsive flares studied bv lSui et ai]d2007| ). Further, 
events were required to have well separated sources to obtain 
spectra of individual sources using imaging spectroscopy with- 
out contamination by other sources. We also analysed partially 
occulted flares since their observation could provide us with 
high-energy resolution spectra of a single source. To achieve that 
we selected partially occulted events with only one source which 
was nearly cospatial at thermal (6-12 keV) and higher (typically 
above 25 keV) energies. Finally, events with a corrected live- 
time below 90% were discarded to avoid pile-up issues. FigureQ] 
shows examples of analysed disc events with several separated 
sources (top) and partially occulted events (bottom). 

Spatially integrated spectra were accumulated during the 
event peak (typically tens of seconds) at energies above 25 keV. 
Spectra of individual sources we re obtained follow i ng the imag- 
ing spectroscopy approach o f iBattaglia & Benz| d2006f) . We 
used the CLEAN algorithm dHurford et al.l 120021) to recon- 
struct images in an interval of several RHESSI spin periods 
covering the time interval for the spatially integrated spectra. 
Using OSPEX0, appropriate regions of interest were chosen 
and feature-integrated spectra, i.e. from a footpoint or a coronal 
source, were obtained. 

3. 1 . Spectral fitting 

Spatially and feature-integrated spectra were forward-fitted 
within OSPEX assuming the following two distinct fitting mod- 
els: 

1 . Isothermal component + thin-target component of an elec- 
tron power-law distribution, i.e. the standard approach. The 
fitted parameters are temperature T and emission measure 
EM of the thermal plasma and parameters of the power-law, 

1 http://hesperia.gsfc.nasa.gov/ssw/packages/spex/doc/ 
ospex_explanation.htm 
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Fig. 1. Examples of events imaged at low (typically 6-12 keV) 
and high (above 25 keV) energies. 

non-thermal, component, i.e. low-energy cutoff E c , spectral 
index S\, normalisation factor nVF, and optionally break en- 
ergy £b and spectral index 62 above the break energy in the 
case of a double power-law. 
2. Thin-target emission of an electron kappa distribution 
(Eq. (f2]i) + two Gaussians to account for iron-line complexes 
at 6.7 and 8 keV. (Note that the contributions of the iron- 
line complexes are included in the isothermal component in 
the standard fitting approach.) The fitted parameters are T K , 
index k, normalisation n^VN K , integrated intensities of two 
Gaussian lines and optionally their centroids E\, Ei. 

A fit was considered acceptable if the value of reduced^ 2 , xl, 
was close to 1 and normalised residuals did not show significant 
clustering at a particular energy. 

The thin-target formula for an electron kappa distribution, 
Eq. (f2j), has been incorporated into the SSW tree and can be ac- 
cessed as OSPEX fitting function thin_kappa . pro . 

4. Consistency of the kappa distribution with flare 
X-ray spectra 

Applying two fitting models described in Section 13.11 to spec- 
tra of several flares, we have found that the spatially integrated 
spectra of disc events cannot be fitted with the electron kappa 
distribution. The best fits using the kappa distribution have large 
value of xl and large residuals, generally around 10 keV. On the 
contrary, those spectra are easily fitted by the standard approach, 
i.e. an isothermal component plus an electron power-law spec- 
trum. Such spectra are mostly a combination of emission com- 
ing from several sources (foopoints, coronal or loop source - 
see Fig. Q] top row) of distinct plasma properties, therefore it 
is difficult to describe the whole spectrum by just one distribu- 
tion. On the other hand, using the combination of an isothermal 
and a power-law component allows us to adjust parameters of 
these functions in such a way that a satisfactory fit usually can 
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14-Nov-2002 11:09:24 - 11:09:52 UT Southern footpoint 



13-Jul-2005 14:13:52 - 14:14:04 UT Coronal source 
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Red: kappa y u = 0.69 
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Hue: kappa y u = 0.67 

n„VN„ T„ k = 0.5357, 0.09, 2.20 
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Fig. 2. Comparison of electron distributions obtained by fitting 
emission from the southern footpoint of the 14-Nov-2002 flare 
(see also Fig. [TJ. Black lines denotes the isothermal plus power- 
law component fit, red and blue lines correspond to kappa distri- 
bution fits. Due to large uncertainties in the data, similarly good 
fits in terms of xl, but differing at energies below ~ 10 keV, can 
be obtained. 



be obtained. The low-energy cutoff helps describe the smooth 
transition between the thermal and non-thermal components in- 
dependently of other fit parameters. Conversely, the shape of the 
kappa function is given only by two parameters: T K and k. 

Next, we examined the feature integrated spectra obtained 
by the imaging spectroscopy approach. Spectra related to foot- 
points can be fitted by both models; the best fits are similar in 
terms of xl an d residuals. The footpoint emission at energies 
above 10 keV is approximated well by the kappa distribution; k 
is generally below 4, therefore the resulting distribution resem- 
bles the power-law used in the standard approach. However, pa- 
rameters T K and n^VN K are not constrained in the model with the 
kappa distribution. Equally good fits can be obtained with values 
of T K and ri^VN K differing by an order of magnitude - compare 
the resulting electron distributions and values of fitted parame- 
ters in Fig. [2] This is probably related to the fact that footpoint 
emission at low energies is reconstructed with high uncertain- 
ties. Therefore, we conclude that there is no advantage in using 
a kappa distribution for the footpoint spectra. 

As concerns coronal sources of disc events, the 13-Jul-2005 
flare (see Fig. [1] top right) coronal emission can be well fitted 
by both spectral models. Both models provide best fits with sim- 
ilar residuals and;^ 2 , close to 1 (see Fig. 01. Moreover, all kappa 
parameters are well constrained. 

Having obtained a reliable fit using the kappa function 
to a coronal source spectrum of a disc event reconstructed 
from CLEAN images, we examined high-energy resolution 
spatially integrated spectra of partially occulted flares. Such 
spectra are generally obtained with much finer energy resolu- 
tion and better signal-to-noise ratio than the data from imag- 
ing spectroscopy (compare uncertainties on the data points in 
Figs. [3]and|4]i. Partially occulted events were recently studied by 
Krucker & Lhi| d2008l) and we analysed several events listed in 
that paper. A representative example is the 30-Apr-2002 event. 
Its spectrum cannot be described well by the model using a 
kappa function (the best fit shows large residuals at energies 
around 10 keV), whereas the standard model matches the spec- 



xl = 1.04 Blue: kappa (n p VN„ T,, k): 0.78, 1.32, 7.96 
Red: Gaussian lines (E,, E 2 ): 6.84, 7.97 
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Fig. 3. Spectrum of the coronal source of the 13-Jul-2005 flare 
(see also Fig. [TJ obtained from CLEAN images. Data are plot- 
ted as crosses; their horizontal sizes correspond to the widths of 
energy bins, vertical sizes represent flux uncertainties. Coloured 
solid lines denote the components of the kappa plus lines fit, 
dashed lines correspond to the components of the isothermal 
plus power-law fit. The thick black line shows the resulting fit 
using the kappa distribution. 



trum reasonably well. This event shows nearly cospatial sources 
at 6 - 12 and 12-25 keV ranges, yet their shapes suggest differ- 
ent structures depending on energy - see Fig.fT] bottom left. 

On the other hand, the spatially integrated spectrum of the 
extended coronal source of 27-Oct-2002 flare seen by RHESSI 
above the limb (see Fig. [1] bottom right) can be fitted well by 
the kappa function up to 50 keV. Again, the standard model pro- 
vides a similarly good fit both in terms of residuals and^J, but 
contrary to the kappa function model a low-energy cutoff must 
be introduced. The best fit using the kappa function is shown in 
Fig.0 

Now the question arises of why the kappa distribution fits 
the spectrum of the October 27, 200 2 coronal source much more 
closely than that of April 30, 2002. iKrucker et al.1 (120071) show 
that the October 27, 2002 source was observed at much higher 
altitudes (0.3 Rq) than the April 30, 2002 event. We think that 
at lower altitudes the X-ray source structure can be more com- 
plex (e.g. several different X-ray sources along the line of sight) 
than for sources at higher heights and thus the fit of the kappa 
distribution is not successful due to the limited number of free 
parameters. 

5. Conclusions 

We have applied a kappa distribution of electrons to RHESSI 
hard X-ray flare spectra. We have shown that spatially integrated 
spectra of disc events, generally composed of contributions from 
coronal and footpoint emissions, cannot be well described as the 
thin-target emission of an electron kappa distribution. On the 
contrary, the standard approach assuming an isothermal and a 
power-law component results in a better fit in terms of residuals 
and reduced if- 2 . This could be due to the fact that those different 
spectral components come from regions of distinct plasma prop- 
erties and therefore it is difficult to describe the spectral shape 
by two parameters, T K and k only. 
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27-Oct-2002 22:53:24 - 22:54:32 UT 



kappa + lines: thermal + power— law: 

xl = °- 77 xl = °' 76 

k = 2.70 6, = 2.14 

E„ = 10.38 




energy [keV] 



Fig. 4. Comparison of the kappa distribution fit (solid) and the 
isothermal plus power-law fit (dashed) to the spatially integrated 
spectrum of the partially occulted flare on 27-Oct-2002 (see also 
Fig. [TJ. Colour coding is the same as in Fig. [3] Top: observed 
photon spectrum (crosses - horizontal and vertical sizes repre- 
sent widths of energy bins and flux uncertainties, respectively) 
and fits (histograms), middle: normalised residuals of the fits, 
bottom: resulting electron distributions. 



Further, from the spectra of individual sources obtained by 
means of imaging spectroscopy, we have found that both foot- 
point and coronal emissions could be fitted using the electron 
kappa function. However, in the case of the footpoints, the fit is 
ambiguous at low energies owing to high uncertainties of data 
in this energy range. Therefore, we cannot conclude whether the 
kappa distribution is appropriate for the footpoint emission. 

In constrast, the results both from imaging spectroscopy and 
high-energy resolution spectra of partially occulted flares sug- 
gest that some coronal sources are consistent with the thin-target 
emission of the electron kappa distribution. 



Summarising the results of all analysed flares, it looks as if 
the decrease of X-ray source complexity increases the relevance 
of spectral fitting by the kappa distribution, especially for coro- 
nal X- ray sources. This is in a greement wi t h the n umerical mod- 
els of lYoon et all (I2006I) and iRhee et all (I2006I) who consider 
particle-wave interactions in plasma and thus describe processes 
appropriate for coronal conditions where Coulomb collisions do 
not dominate. 

Our results do not rule out the widely used interpretation of 
hard X-ray spectra being produced by thermal and power-law 
components, yet they show that other distributions, such as the 
kappa distribution, are consistent with the RHESSI data as well. 
In our view, the advantage of using the kappa distribution lies 
mainly in the facts that the kappa distribution is a natural con- 
sequence of physical processes in space plasma, that such a dis- 
tribution is applicable to in-situ detected electrons, and no low- 
energy cutoff in the mean electron spectrum has to be consid- 
ered. 

We note that the low-energy part of hard X-ray spectra is 
important to constrain the best fit. For future studies we pro- 
pose to combine hard X-ray data with a detailed analysis of 
soft X-r ay continuum and line s from other instruments such as 
RESIK dSvlwester et al.ll2005l) or the recently la unched SphinX 
on board Coronas-Photon dSvlwester et al.ll2008l). Fina lly, on the 
basis of nonextensive thermostatistics. lLeubnerl d2004) proposed 
a so-called double-kappa distribution and applied it to observed 
interplanetary electron distributions. Future RHESSI data anal- 
yses may use this idea as an alternative approach to the single 
kappa distribution used here. 
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